Time-of-arrival distribution for arbitrary potentials and Wigner's time-energy 
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A realization of the concept of "crossing state" invoked, but not implemented, by Wigner, allows to 
advance in two important aspects of the formalization of the time of arrival in quantum mechanics: 
(i) For free motion, we find that the limitations described by Aharonov et al. in Phys. Rev. A 
57, 4130 (1998) for the time-of-arrival uncertainty at low energies for certain measurement models 
are in fact already present in the intrinsic time-of-arrival distribution of Kijowski; (ii) We have also 
found a covariant generalization of this distribution for arbitrary potentials and positions. 
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In spite of the emphasis of quantum theory on the con- 
cept of "observable", the formalization of "time observ- 
ablcs" is still a major open and challenging question. The 
"arrival time" has in particular received much attention 
in the last few years (see for a recent review) . Consid- 
ering several candidates proposed for the time-of-arrival 
distribution in the simple free motion, one dimensional 
case, some of us have recently argued || in favor of a 
distribution originally proposed by Kijowski 11^, be- 
cause it satisfies a number of desirable conditions. 

This distribution can be associated with a POVM and 
obtained in terms of the eigcnfunctions |T, a; X) (a = ±) 
of the time-of-arrival operator T introduced by Aharonov 
and Bohm 0] 

f=~[(q-X)l+i(q-X)}, (1) 
2 P P 

U K (t + T; X- = E K T ' a = X \^))\ 2 . ( 2 ) 

a 

(We consider here the general ID case with both positive 
and negative momenta as in [^]|^||.) where m is the 
mass, X is the arrival point, q and p are position and 
momentum operators, and 



\T, a; X) = e maT / h {\p\/m) l / 2 Q{ap)\X) . 



(3) 



(The operator |pp/ 2 is defined by its action on mo- 
mentum plane waves, IpI 1 ^ 2 ^} = IpI 1 / 2 ^}.) Urt^ + 
T; X;ip(t)) represents the probability density of arriv- 
ing at X, at the instant T + 1, for a given wavepacket 
if) (t), and t is the parametric time that characterizes the 
evolution of the state tp(t). (Typically one sets t = so 
that T is the "nominal arrival time".) This distribution 
satisfies in particular the important covariance condition 
under time translations, n^-(i + T — t'\ X; ijj(t + t')) = 
Tl K (t + T; X; ip(t)). For other properties see 

In this letter we study two major aspects of this distri- 
bution that had not been addressed: (z) For states with 



positive or negative momenta we shall obtain the states of 
minimum time uncertainty (for given energy width), and 
find the same type of limitation pointed out by Aharonov 
et al. ||; (ii) we shall also generalize (||) for arbitrary po- 
tentials. 

To handle conveniently this two issues let us first elab- 
orate on the form of Hk- For a = + the contribution in 
(||) can be interpreted as a quantum version of the posi- 
tive flux at the time T + 1 due to right moving particles. 
Similarly, for a = — , one has a quantum version of mi- 
nus the negative flux due to left moving particles, again 
a positive quantity. Explicitly, 

Tl K , a (t + T) = ty(t + T)\(\p\/m)^ 2 Q(ap)6(q- X) 

x e{ap)(\p\/m) 1/2 \^(t + T)} , (4) 

where 5(q — X) = \X)(X\. The positive operator in (^) 
corresponds to the classical dynamical variable 

S(q-X)^Q(ap), (5) 

TO 

whose average represents the modulus of the flux of par- 
ticles of the classical ensemble that arrive at X from, one 
side at a given time. This connexion was pointed out by 
Delgado |@. There are of course many possible quan- 
tizations of this quantity but the symmetrical one in (^) 
turns out to be the only one that satisfies the symmetry 
and minimum variance properties of 11^. 

It is useful to write the positive operator in ([|) in the 
form \u a )(u a \, where 



\u a ) = \T = 0,a;X) 



(6) 



is the "crossing state". As emphasized in it is im- 
portant to keep in mind that, being non-normalizable, 
its literal interpretation is problematic. Only normalized 
wave packets peaked around these states have proper- 
ties as close as desired to the sharp crossing behaviour 
expected on intuitive grounds. 
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Let us first discuss the point (i) related to the time- 
energy uncertainty principle. Since the Hamiltonian and 
the time operator (Q) are conjugate variables a minimum 
uncertainty product can be established in the usual fash- 
ion However, Aharonov et al. have proposed, based 
on a series of models, a second limitation on the possi- 
ble values of the time-of-arrival uncertainty: St Ek > h~, 
where St is the width of the "pointer variable" used to 
measure the arrival time, and Ek "the typical initial ki- 
netic energy of the particle" j|] . It is to be stressed that 
this relation is based on measurement models for the ar- 
rival time where some extra (clock) degree of freedom is 
coupled continuously with the particle. We shall see that 
in fact the "intrinsic" distribution Tlx (there is no explicit 
recourse to additional pointer degrees of freedom to de- 
fine TIk) is consistent with the behaviour that Aharonov 
et al. described for their models 

There are many other time-energy uncertainty rela- 
tions but here we shall be mainly interested in the 
version of E. P. Wigner Jl2| , because his formalism is par- 
ticularly suited for the time of arrival. In the original pa- 
per Wigner did not consider in detail any application and 
described a variational method to find the states of min- 
imum uncertainty product, but did not actually obtain 
these states, except in two analytically solvable cases. 
We shall extend Wigner's work in several directions by 
evaluating the states of minimum uncertainty product 
and applying the formalism to the arrival time. For com- 
pleteness we shall next briefly summarize the main re- 
sults obtained by Wigner in ]l^ |, and add a number of 
comments and observations relevant for our application. 

He defines the basic amplitude as \(t) = (H^W)' 
where \u) is in principle any state vector. (Wigner's 
formalism encompasses many different time-energy un- 
certainty products, depending on the \u) chosen, each 
with its own physical interpretation.) Note that t is 
considered the independent variable of x, and \u) is 
fixed. \u) is not necessarily a Hilbert space normaliz- 
able vector. It may also be, for example, a position, 
or a time-of-arrival eigenvector provided the minimum 
uncertainty product state obtained is square integrable. 
P(t) = \x{t)\ 2 1 \x{t)\ 2 dt plays the role of a normal- 
ized distribution for being found at \u) at time t. This 
is not a standard quantity in the ordinary formulation of 
quantum mechanics (which assigns probabilities only at 
fixed time i), but the interpretation is consistent with the 
ordinary formulation in the following way: Here "being 
found" implies operationally to perform measurements of 
at a given time t for the members of an ensemble 
prepared at to < t according to ip^o)- This is repeated at 
different times but the ensemble is always prepared anew 
at to with the same specifications. The distribution of 
positive counts as a function of time is proportional to 
|(u|-0(i))| 2 , and P(t) is obtained when this distribution 
can be normalized (which is not always the case). (P{t) 
does not correspond to a continuous measurement that 



modifies tp(t)-) The moments of P(t) are defined in the 
usual way, and in particular the second moment with re- 
spect to to is defined as 

/oo />oo 
\ X (t)\ 2 (t-t a fdt/ \ X (t)\ 2 dt. (7) 
-oo J —oo 

The information contained in %(£) can also be encoded 
in its Fourier transform, 

r)(E) = h- 1 ' 2 J x(t)e lEt/n dt. (8) 

The conjugate variable E has dimensions of energy, but 
rj{E) is not, in general, an energy amplitude in the con- 
ventional sense. r\ and the conventional energy amplitude 
of ip(t — 0), can be related by expanding tp(t) in a basis 
of energy eigenstates \E, a), 

r)(E) = h 1 ' 2 (u\E, a) (E, a \^{Q))@{E) . (9) 

Oi 

In a general case a is an index to account for the possible 
degeneracy. In particular, for free motion, a = ±, and 

lE ' a) = \2E) \P = a ( 2mE ) 1/2 ) ■ ( 10 ) 

Analogously to ((?]) the second energy moment with re- 
spect to Eq is defined as 

/•oo />oo 

e 2 ^ \ V (E)\ 2 (E-E ) 2 dE/ \ V (E)\ 2 dE. (11) 
Jo Jo 

Neither to nor Eq should in general be identified with the 
average values of time and energy. They are instead ref- 
erence parameters fixed beforehand to evaluate the mo- 
ments. As a consequence, e 2 and r 2 should not in gen- 
eral be identified with the "variances" (AE) 2 and (At) 2 , 
which are the second moments with respect to the aver- 
age values. 

Since rj(E) and x(t) are Fourier transforms of each 
other the uncertainty product er is greater than fi/2 (a 
peculiarity of time and energy with respect to position 
and momentum is that the equality is never satisfied 
p2|). In fact the bound increases substantially as Eq 
decreases. Wigner sought for the function rj(E) that ren- 
ders r to a minimum for fixed e. In order to have a 
finite second moment r 2 , r](E) must vanish at the origin, 
77(0) = 0, so 77 must vanish at both ends of integration 
E = 0, 00. Using partial integration, and the notation 
770 = rier lEt °l h (Wigner showed that the minimum of re 
must correspond to a real 770.) , one finds that 

/•oo />oo 

r 2 = h 2 \dT] (E)/dE\ 2 dE/ \ Vo {E)\ 2 dE. (12) 
Jo Jo 

The product r 2 e 2 subject to the constraint of fixed e 2 is 
minimized by variational calculus. This leads, using ( |Tl| ) 
and @, to 
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£o) 2 r?o - (r 2 + A')?7o , (13) 



where A' is a Lagrange multiplier. This equation is for- 
mally similar to the Schrodinger equation for the har- 
monic oscillator, except for the boundary condition at 
E = 0, 770 (0) = 0, and the subsidiary condition for e (11). 
The minimum r is obtained from the lowest eigenvalue 
corresponding to the value of A' where the subsidiary con- 
dition is satisfied. Fortunately the solution depends only 
on the ratio E /e, namely i](E;E ;e) = g(E /e; E /e), 
where g is the solution of ( |l3| ) with e = 1. Note also 
that, since |?7o| 2 = r/| 2 , the value of to does not play any 
role in the minimization process. (Physically the state 
Xo(t) corresponding to to = is valid for any other time 
to by a shift of the argument.) The minimization of r 
for fixed Eo/e requires a method to solve the differential 
equation ( Jl3| ) for many different values of A', until e 2 = 1 
is satisfied. In our calculation the successive values A' (n) 
have been obtained with a Newton-Raphson algorithm, 
and the lowest eigenstate and eigenvalue of (|l3| ) for each 
A' (n) with a very efficient "relaxation method" |l3| . 

The only explicit case considered in the original paper 
by Wigner was \u) = \X). The corresponding P(t) pro- 
vides a time distribution for the presence of the particle 
at X but not for its arrival. In an intriguing "general 
observation", Wigner stated that, instead of asking for 
the probability that the particle be at a definite land- 
mark in space, just at the time t, "it would be more 
natural to ask ... for the probability that the particle 
crosses the aforementioned landmark at the time t from 
the left, and also that it crosses the landmark, at a given 
time, from the right." But the paragraph ends with "This 
point interesting though it may be, will not be elab- 
orated further" . Precisely, it is our aim here to elabo- 
rate further this question. Indeed the probabilities men- 
tioned by Wigner can be obtained by means of the cross- 
ing states \u+) and discussed before, see (||). Spe- 
cializing to states having only positive/negative momen- 
tum Hk = W.K,a(t) — \{u a \%)j(t))\ 2 provides the arrival 
time distribution at t. Considering %(t) = (u a \ip(t)), 
we see that Wigner's probability density is nothing but 
Kijowski's distribution, P(t) = U K ^ a (t) = |x(i)| 2 . More- 
over, the Fourier transform of x(t) is i n this case the 
standard energy amplitude, f](E) = (E, a\ip(0)), so that 
e becomes the spread (around Eq) of the ordinary energy 
distribution. 
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FIG. 1. er (in atomic units, h — 1) versus {E)/e for the 
states of minimized uncertainty product (dashed line); e/{E) 
(solid line). 

For a given ratio Eo/e there is a minimum value of re. 
The family of states of minimal uncertainty r](E;Eo;e) 
with the same ratio Eq / e have in common the same value 
of (E)/e. Figure I shows re versus (E)/e for the states 
of minimized time energy uncertainty product. For com- 
parison we also show the curve e/{E). Clearly 



t > h/(E) , 



(14) 



which has the same form as the relation proposed by 
Aharonov et al. based on measurement models 0. How- 
ever, r is not due to the effect of any measuring appa- 
ratus, it is an intrinsic uncertainty associated with an 
intrinsic time-of-arrival distribution. It is not the cou- 
pling introduced in a measurement between the particle 
and other degrees of freedom that leads to this relation 
but the very quantum mechanical nature of the particle 
alone and the lower bound of the energy. 

To elaborate the figure, Eo/e has been increased regu- 
larly from the minimum possible value —1. (For smaller 
values it is impossible to satisfy the subsidiary condition.) 
For each value the minimization of re is performed and 
the corresponding (E)/e is obtained. As (E)/e — + 00 the 
minimum uncertainty product tends to the (global) mini- 
mum h/2, the same value found for position and momen- 
tum, because the effect of the lower bound of the energy 
tends to disappear in that limit, and r)(E) becomes closer 
and closer to a Gaussian centered at Eq with variance e 2 
However, in the opposite limit the lower bound at 
E = plays an important role. Indeed, (E)/e — > corre- 
sponds to the limit Eo/e — > —1, and the only way to sat- 
isfy the constraint is by strictly localizing rj(E) at E = 0, 
(AE — > 0), but this completely delocalizes the conjugate 
time variable, namely r — > 00. Thus, (14) appears as a 
consequence of the ordinary uncertainty principle, due to 
the tendency of the minimum uncertainty product states 
to have smaller variances for smaller energies. (The pre- 
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cise dependence for arbitrary values of Eq/c has to be 
obtained numerically.) 

The second question we shall address is the general- 
ization of the free motion distribution (^) for arbitrary 
potentials and positions. A generalization based on a 
quantization of classical expressions as in (|l]) is problem- 
atic: the classical expressions for the time of arrival will 
rarely be analytical, not all phase space points lead to 
arrival, and the ordering problems may be formidable. 
The way out though, will be surprisingly simple in terms 
of crossing states. 

There is in fact nothing that limits (^]), and the cor- 
responding operator in (|4|) to free motion. In particu- 
lar, the classical motivation for considering \u a ) a "cross- 
ing state" is equally valid when an arbitrary potential is 
present. The state (||) may be regarded as one that has 
evolved "backwards" with Ho a time T from the crossing 
state \u a ) so that it becomes \u a ) precisely at the nomi- 
nal arrival time T . In the same vein we construct for an 
arbitrary Hamiltonian H 



It is evidently covariant under time translations as it 
should; in general it is not normalized, and it may be un- 
normalizable (its classical counterpart shares these prop- 
erties). For example, it may be constant for stationary 
states, or periodic for oscillating coherent states in a har- 
monic potential, but it provides in any case relative in- 
formation by comparison of two times. Consistent with 
this, the states \T, a;X) do not form in general a com- 
plete basis. According to its classical analog it takes into 
account any crossings (not only first, or last). To illus- 
trate the distribution we have evaluated II + for a colli- 
sion with tunnelling at three different positions before, in, 
and after a square barrier, see Figure 2. The only previ- 
ous attempt to generalize Hk applied only to asymptotic 
positions where the motion is essentially free fll4| . Our 
generalization is instead valid for arbitrary positions. 

One of us (J.G.M.) acknowledges C. R. Leavens for 
helpful discussions. This work was supported by Go- 
bierno Autonomo de Canarias (PB/95), MEC (PB97- 
1482), and CERION. 



\T, a ;X) = e iHT ' n \u a ) 



(15) 



so that at the nominal arrival time T, \u a ) is recovered. 
The generalization of the arrival time distribution for ar- 
bitrary potentials is therefore 



II(t + T; X; </>(*)) = J2 \(u a \e~^ T ' h m)W 



(16) 



n 




FIG. 2. Current density and 11+ (t) before the barrier 
(X — 8, solid line and circles), at the barrier (X = 12.25, 
dashed line and squares), and after the barrier (X = 15, long 
dashed line -in this case the two functions are indistinguish- 
able in the scale of the figure-), for a collision of a Gaussian 
state (with minimum position-momentum uncertainty prod- 
uct at t = and (q) = (p) = 5, ((q — 5) 2 ) = 1) of a particle 
with mass m — 1/2 versus a square barrier of energy 40 lo- 
cated between x = 12 and 12.5. All quantities in atomic 
units. 
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